Stability of curvature perturbation with new covariant form 
for energy-momentum transfer in dark sector 



Cheng- Yi Sun*' 1 , Yu Song 1 and Rui-Hong Yue 2 

institute of Modern Physics, Northwest University, Xian 710069, P.R. China. 
2 

Faculty of Science, Ningbo University, Ningbo 315211, P.R. China. 

October 11, 2011 



Abstract 

It was found that the model with interaction between cold dark matter (CDM) and dark 
energy (DE) proportional to the energy density of CDM p m and constant equation of state of 
DE Wd suffered from instabilities of the density perturbations on the supper-Hubble scales. 
Here we suggest a new covariant model for the energy-momentum transfer between CDM 
and DE. Then using the covariant model, we analyze the evolution of density perturbations 
on the supper-Hubble scale. We find that the instabilities can be avoided in the model with 
constant Wd and interaction proportional to p m . Furthermore, we analyze the dominant 
non-adiabatic mode in the radiation era and find that the mode grows regularly. 

1 Introduction 

We are convinced by the increasing observations [H El E] that the present universe is dominated 
by the so called dark energy (DE), which accounts for ~ 70% of the critical mass density and has 
been pushing the universe into accelerated expansion [H [5] . And the other main component in 
the universe is cold dark matter (CDM), which accounts for ~ 30% of the critical mass density 
and behaves as the pressureless dust. Then it is natural for us to consider that the two dark 
components might interact mutually. And furthermore it is found that an appropriate interaction 
can help to alleviate the coincidence problem [6j [7] , namely why DE and CDM are comparable 
in size exactly today |5J. Different interacting models of dark energy have been investigated 
intensively [8j |9] . 
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Usually, in the literature, the model with interaction between DE and CDM is described by 
the two continuity equations 

p m + 3Hp m = Q, (1) 
p d + 3H(1 + w d )p d = -Q, (2) 

where Q denotes the phenomenological interaction term between DE and CDM; p m and p d are 
the energy densities of CDM and DE respectively; w d = p d /p d is the equation of state parameter 
of DE; p d is the pressure density of DE; H = a/a is the Hubble parameter; a(t) is the scale factor 
in the Friedmann-Robertson- Walker (FRW) metric; a dot denotes the derivative with respect to 
the cosmic time t. In the note we do not allow the phantom case w d < —1. Owing to the lack of 
the knowledge of micro-origin of the interaction, usually the interaction term is parameterized in 
a simple form as [6] 

Q = 3H(ap m + (3p d ), (3) 

where a and (3 are positive constants. The interaction term Q would influence not only the 
background dynamics of the universe, but also the growth of the perturbations of the cosmological 
fluids. 

Recently, in Ref.[10j, by modeling DE as a fluid with constant w d > —1, the authors investi- 
gated the evolution of the linear density perturbations and have shown that the combination of 
constant w d and the simple interaction form Q given in Eq.(|3]) leads to an instability: the cur- 
vature perturbation on the super-Hubble scales blows up in the early universe jlOJ . The explicit 
models investigated in Ref. [10J included the two cases of = and a = (3 in Eq.(j3]). Further 
more, in [12], it is concluded that the perturbations in the dark energy become unstable for any 
model with constant w d > — 1 and non-zero a, no matter how small the parameter a is made. In 
PHI [121 IE], the case of a = was surveyed and it was found that the instability can be avoided 
if is made small enough. In [14] . by modeling DE as a quintessence field, the author found that 
the instability can also be avoided even for the interaction Q proportional to p m . 

Then it seems that the model with constant w d and non-zero a in Eq.([3]) is ruled out as a 
viable interacting model. However, in the note we try to show that the evolution of the density 
perturbations becomes regular even in the case of /3 = in Eq.([3]) if we adopt a new covariant 
form for energy-momentum transfer between DE and CDM. 

The remain part of the note is organized as follows. Firstly, we display our new covariant 
model for interaction in the dark sector. Secondly, by assuming the universe filled only by DE 
and CDM, we investigate the evolution of the density perturbations and show that instability 
can be avoided. Thirdly, by considering the effects of the radiation (photons and neutrinos) 
and baryons, we survey the dominant non-adiabatic mode in the radiation era and show that no 
non-adiabatic mode blows up. Finally, conclusions and discussions are given. 
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2 A Covariant Model For Dark- Sector Interaction 



The conservation laws ([I]) and (|2J) work well in describing the background evolution of the universe. 
But in order to study the evolution of the density perturbation, we need a covariant form for the 
energy-momentum transfer between DE and CDM which holds in an inhomogeneous universe 
and reduces to Eqs.flT]) and (j2J) in a homogeneous FRW universe. Usually, the covariant form for 
energy-momentum transfer is taken to be [T5l HD] 

v*ir=<& E^ = o, (4) 

A 

where A — m, d to denote CDM and DE respectively. In Ref.[10j, it is assumed that 

Qm = ~Qd = Q<, (5) 



where is the four velocity of CDM. The conclusions in [lOj are based on the above ansatz. 

However, we find that there may exist the other natural covariant form for the energy- 
momentum transfer. Let us show it. In the note, we take the interaction Q to be 

Q = 3aHp m , (6) 

This is just the case of = in Eq.Q. Then we can rewrite Eqs.(pQ) and (j2J) as 

p m + 3H{p m + p e £) = (7) 

p d + 3H(p d + pf) = 0, (8) 

where 

Pm = -ap m , (9) 

pf = p d + ap m . (10) 

Motivated by these equations, we may define the effective energy-momentum tensors of CDM and 
DE respectively as 

■Ci = PmU^ n u m + p m {u^ n u m + ), (11) 

= Pd« + PfW + O. (12) 

where ut^ and are the four velocities of CDM and DE respectively. And the two effective 
energy-momentum tensors are conserved respectively 

TZ-„ = T ed;tM = (13) 
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It can be easily checked that Eqs.flT]) and (j2j) can be reduced from Eg. (TT3l) in the background FRW 
universe. The corresponding Einstein equations should be 

G^ = 8irG(TZ + T£). (14) 

Now we can survey the evolution of the curvature perturbation by expanding Eqs. (113j) and (TT41) 
to the first order of the density perturbations. 

The covariant model given in the last paragraph is very different from the one defined in 
Eq.flU), although both of them give the same background evolution of the universe. In the model 
in Eq.pj, the perturbations of the interaction effect the evolution of the density perturbations 
via the continuity equations and do not appear in the Einstein equations explicitly, while in the 
model defined in the last paragraph, the perturbations of the interaction appears explicitly both 
in the continuity equations and in the Einstein equations. 

3 Evolution of Density Perturbations 

In the section, we apply the covariant model defined in Eqs.( |T3i) and ( |T4"1) to study the evolution 
of the density perturbations in the model with Q given in Eq.([6]) and constant Wd- For simplicity, 
we consider a flat FRW universe filled only by DE and CDM. The perturbed FRW metric in the 
conformal Newtonian gauge is given by 

ds 2 = a 2 (r)[-(l + 2<p)dr 2 + (1 - 2ip)dx 2 }, (15) 

where <fi and ip denote the scalar perturbations. The corresponding Friedmann equation can be 
rewritten as 

o {a'\ 2 8ttG , . 9 

W = -^-^ + Pd) a - ( 16 ) 

Hereafter, primes denote the derivatives with respect to the conformal time r. With Eq.([6]) and 
constant Wd, the background continuity equations ([T]) and ([2]) can be solved exactly: 

p m = p m0 a- 3{1 - a) , (17) 
Pd = p d oa-^ 1+Wd) + (^-) Pm oa-\a-^ - a 3a ). (18) 

Hereafter, the subscript denotes the present value of the corresponding parameter and ao = 1- 
3.1 Evolving Equations of Density Perturbations 

When the perturbed metric in Eq. (jl5l) is considered, the four velocities of CDM and DE are 

< = a~ 1 (l-0,«9 i t; m ), = 0-^1 -0,9^), (19) 
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where v m and v d are the peculiar velocity potentials of CDM and DE respectively. Usually, we 
define the volume expansion rates of CDM and DE (in Fourier space) respectively as 

m = -k 2 v m , 9 d = -k 2 v d . (20) 

We use Sp m , Sp d , Sp d , 5p e ^ and 5p c d to denote the first-order perturbations of the corresponding 
parameters. Then we can introduce the dimensionless factional density perturbations of CDM 
and DE as 

x &Pm r 5p d 

m = , d = • {21) 

Pm Pd 

The curvature perturbation on the constant-p^ (A = m, d) surface and the total curvature per- 
turbation on the constant-p (p = p m + pd) surface are defined respectively as 

> , njbpA y. , aj Sp m + 5p d , . 

Pa Pm + Pd 

From Eq.()n]), we have 

6p* = -ap m 5 m . (23) 

And from Eq.( lTU|) . we have 

$pf = 5 Pd + ap m 5 m . (24) 
Here, following the analysis in Ref.[10j, we take 

S Pd = 5p d + (1 - w d )[3U{\ + w d + a— )p d ]§. (25) 

Pd K 

Then from Eq. fU3l) . we can get the evolving equations of the density and velocity perturbations 
of CDM and DE as 

6' m + (1 -a)9 m - 3(1 - a)ij' = (26) 

6' m + H(l + 3a)6 m -k 2 <J) + k 2 -?—5 m = (27) 

I — a 



and 



8' d + 3H(1 -w d - a^)5 d + 9H 2 (1 - w d )(l + w d + a^)% 
Pd Pd k 2 



+3aH—5 m + (1 + w d + a— )(8 d - 3^') = 
Pd Pd 



(2? 



a(l — a)— a— A 
9' - 2W d - 3H— -^j-6 d - k 2 (f) - k 2 e ^^S m - k 2 - — = (29) 
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In the conformal Newtonian gauge, the first-order perturbed Einstein equations Eq. (fT4l) give 
us pi] 

3W + k 2 ip + 3V. 2 (f) = -AnGa 2 (5 m p m + 8 dPd ) , (30) 
&V + k 2 U<P = AnGa 2 [(1 - a)p m 6 rn + (1 + w d + a^)p d 9 d ] (31) 

Pd 

r + %W + <P') + (2- - H 2 )^ + ^(V - <f>) = 4vrGaW m ff + dpf ) (32) 

a 3 

^-0 = 0. (33) 

Only two of the above equations are independent. Choosing two of them (e.g. Eqs.( l30l and 
(|33|) ). and using Eqs. (l26|) . (1271) . (128]) and (|29|) together, we can solve these evolving equations 
numerically, if the initial conditions are given. 



3.2 Adiabatic Initial Conditions 

In the early universe, a C 1, Eqs. f|T7|) and ffTgj) indicate 



en -> -!^±^, (34) 



and then, from Eq. (fl6l) . we have 



U = T ^-T~\ r= - l— JY!*+!La\W (35) 



1 - 3a (1 - 3a)H \ w d il m0 

where 

"mO — • 

Here we adopt the adiabatic initial conditions to study the evolution of the density perturbations 
on the super-Hubble scales (k <C aH). To the lowest order in kr, we can set the adiabatic 
conditions as 

(f> = ip = = Const., (36) 

8m = $d = (37) 

w d + a 

(w d + a + aw d - a 2 )(l - 3a) 2 
9m = 3(1 -a 2 )(w d + a) k (38) 

d 3(1 -a)(w d + a) (2w d + a -1) V ; 
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Fig. 1: log 10 |C| versus log 10 a in the interacting model for fixed fl m o = 0.3, h = 0.67, k = 
1.5 x 10~ 4 Mpc _1 , a = 10~ 3 and different w d . 

3.3 Evolution of Curvature Perturbation 

Now using the initial conditions given in the last subsection, we can solve the evolving equations 
fl26|) - fl30|) and fl33|) numerically to obtain the evolution of the density perturbation and then get the 
evolution of the curvature perturbation £ by using Eq.f l22|) . We have displayed the results in FigJT] 
and FigfSJ In the two figures, we have fixed the parameters as f2 m0 = 0.3, k — 1.5 x 10~ 4 Mpc -1 , 
Hq = 100/i km sec _1 Mpc _1 and h = 0.67. And we have taken = 10~ 25 and the initial values 
of a to be a = 10 -11 . In FigHJ we show the evolutions of log 10 \(\ for fixed a = 1 x 10~ 3 and 
different Wd- And in Figj2]we show the evolutions of log 10 \(\ for fixed Wd = —0.94 and different 
a. The evolutions displayed in FigJT] and Figj2] manifest the standard power-law growth and no 
instabilities are present, which are similar the results in |T4"] . 

4 Dominant Non-Adiabatic Mode 

In the last section, we have shown that in our new covariant model, instabilities on the super- 
Hubble scale can be avoided. But the conclusion is obtained by assuming that the universe is 
filled only by CDM and DE. In this section, we discuss the dominant non-adiabatic mode deep in 
the radiation-dominated era by including the components of photons, neutrinos and baryons in 
the universe. If the dominant non-adiabatic mode evolves regular in the radiation era, we believe 
the instabilities can also be avoided even when the components of radiations and baryons are 
involved. 
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Fig. 2: log 10 |C| versus log 10 a in the interacting model for fixed fl m0 = 0.3, h = 0.67, k = 
1.5 x 10 -4 Mpc -1 , w d = -0.94 and different a. 



The v4-fluid energy-momentum tensor including perturbations is taken to be 

Tf = (Pa + Va)u\u\ + PAS T + <, (40) 

where u\ is the four velocity 



a 1 (l - (f),diV A y 



We have allowed an anisotropic shear perturbation ir 1 ^ , and A = b, 7, v to denote the correspond- 
ing parameters of baryons, photons and neutrinos. We take tt^ u = ir!j u = and 



7r°f = 0, Trj? = a" 2 ($0,- - -<^>,. (41) 

The fluids of CDM and DE are described by the effective energy-momentum tensors defined in 
Eqs.f lTT]) and ffT2]) respectively. 

Early in the radiation era, the Friedmann equation reads 

H 2 a- 2 = —(p 7 + P„) = ^-Proa" 4 - (42) 



Then we have 

^tfor, H^r- 1 , n r0 ^^^. (43) 

3ii 
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In the radiation era, the perturbed Einstein equations give us that 

3r -y + k 2 ip + 3r" 2 = -4nGa 2 (s Pm + 5p d + J2 <W) , (44) 

A 

k 2 (ip' + r"V) = ^Ga 2 [(p m + p^)9 A + (p d + pf)9 A + ^( P a + Pa)#a] , (45) 

A 

r + 2r-y + r-y - r" V + ^(V - 0) = 4vrGa 2 (V m ff + 5pf + V) 5p A ) , (46) 
^ - = SttGtt^, (47) 

where A runs over 6, 7, and za 

For CDM and DE, the perturbed continuity equations are given by Eqs. (!26|) - (!29|) . For baryons, 
the perturbed continuity equations (in Fourier space) are [ID] 



and for photons [ID] 



and for neutrinos ITD 



5' b = -9 b + 3^', (48) 
9' b = -W b + k 2 <f), (49) 



<j; = ~0 7 + 4^, (50) 
0' = \k 2 5 y + k 2 (J), (51) 



5' V = ~9 V + ^, (52) 

^ = + k 2 <p - k 2 a u , (53) 

< = ^0 U , (54) 

where 9 A = —k 2 VA and 0^ = 2k 2 7i„/[3a 2 (p u + p„)] 

In order to find the dominant non-adiabatic mode, we assume a leading-order power law for 
perturbations [TDJ 

V = ^(fcr)"*, = Mkr) n *, S A = B A {kr) n \ 9 A = C A (kr) SA , a v = D v (k,r) n ° . (55) 
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Here the subscript A — c, d, b, 7, and v denotes the corresponding parameter of CDM, DE, baryons, 
photons and neutrinos respectively. To the leading order in kr, the equations ()26l) -(l29l) and ( 14^1) - 
fl54l) may be solved, in terms of ip: 

= J ^' J = 1 -5K + 2)Mi ) + 8ij; (56) 

<L = S v = 6„ = 9 V = ^\k 2 ri), (57) 

7ty + 1 

^ = 3^, ^ = _^—A; 2 r^, (58) 

5 m = 3(1 - # m = ~ a _ fcVV, (59) 

+ 2 + oa 

5 d = r0 1+3a (w d + a)K + J + 2)-^, (60) 
a\L m0 n t 

9(1 - w d ){l - a) 



where R v = p v j (p 7 + p y ) and 



n t = j^BHl, (62) 

Eq.( )60|) indicates that the modes are regular (i.e. well behaved as kr — > 0) provided 

Re[n^] > 1 + 3a. 

For Wd ~ — 1, this leads to 

3 1 
— Wd > 1 + 3a =>- a < ■—. 
2 ~ 6 

Correspondingly, the total curvature perturbation ( is defined as 

22 a Pa 

where A runs over m, d, 7, v and b. Then £ can be expressed in terms of ip as 

C = --K + J + 2)^. (63) 
For u><j ~ — 1, n,/, is a complex number and 

r 1 3 
Re[r^] ~ -. 

So the dominant non-adiabatic mode grows in a regular power law and no instabilities are present. 



10 



5 Conclusions and Discussions 



In the note, we have suggested a new covariant model for dark-sector interaction to avoid the 
instabilities of the curvature perturbation on the supper-Hubble scale. By using the covariant 
model and choosing Q = 3aHp m and constant w^, we analyze the evolution of density perturba- 
tions in a universe filled only by CDM and DE and find that the instabilities shown in pU] can 
be avoided in our covariant model. Further more, we analyze the dominant non-adiabatic mode 
early in the radiation era, and find the non-adiabatic mode evolves regularly. So we believe that in 
our covariant model, the instabilities can also be avoided even in the universe filled by radiation, 
matter and DE. 

Actually, it is not very surprising for the instabilities to be avoided in our covariant model. We 
know, from Eqs.(TTTT)-(fHl). that the two mutually-interacting dark fluids can be described respec- 
tively by two effective energy-momentum tensors which are conserved separately. This means the 
interacting model can be taken as a non-interacting model effectively, while no instabilities are 
present in non-interacting models. So we expect the instabilities can be avoided in our covariant 
model. 
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